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Figure  1.1 


A  Typical  Plot  of  Stability  Regions 


PLAN  OF  THE  REPORT 


There  are  four  sections  of  this  report. 

1.0  This  gives  a  general  overview  of  what  Is  known  and  what  Is 

commonly  done  for  the  general  viscous  problem. 

2.0  To  obtain  a  qualitative  feel  for  what  may  be  achieved  we  turn 

In  section  2.0  to  the  Inviscid  problem.  The  general  statements  that  can 
be  made  about  stability  when  boundaries  are  modified  are  described. 

3.0  This  section  is  devoted  to  specific  inviscid  examples.  It  is 

shown  that  dramatic  changes  in  stability  are  achievable. 

4.0  In  section  4.0  an  interesting  stability  problem  which  arises  in 

boundary  layer  theory  is  analyzed. 


vii 


1.0 


GENERAL  DISCUSSION 


1 . 1  Introduction 

The  drag  force  on  a  body  moving  through  the  water  is  largely 
determined  by  the  nature  of  a  thin  boundary  layer  between  the  stagnant 
fluid  at  the  body  surface  and  the  free-stream  fluid.  At  low  speeds  the 
flow  in  the  boundary  layer  is  smooth  (laminar)  and  the  drag  is  low. 
Beyond  a  well-defined  critical  speed,  the  flow  becomes  chaotic 
(turbulent)  and  the  drag  is  much  larger. 

To  reduce  drag,  one  must  manipulate  the  boundary  layer  in  such 
a  way  that  either  the  critical  speed  for  transition  to  turbulence  is 
increased  beyond  the  desired  operating  speed  or,  given  that  the  flow  is 
turbulent,  the  turbulent  flow  is  modified  in  a  way  which  reduces  the 
actual  drag  force.  The  classic  motivation  to  study  this  problem  is  of 
course  to  Increase  the  maximum  speed  of  a  vehicle  with  a  given  power 
plant.  A  second  motivation,  relevant  to  SSBN  security,  is  the  desire  to 
reduce  the  intensity  of  the  water  disturbance  caused  by  the  passage  of 
even  a  slowly  moving  vehicle. 

Since  turbulence  is  an  extremely  complicated  and  ill-understood 
phenomenon,  we  will  not  consider  the  problem  of  manipulating  the 
turbulent  boundary  layer,  despite  its  overwhelming  practical 
importance.  Our  remarks  are  entirely  directed  to  the  problem  of 
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Increasing  Che  maximum  speed  aC  which  the  boundary  layer  can  remain  in 
laminar  flow,  a  problem  which  can  be  studied  using  the  well-explored 
mathematics  of  linear  stability  theory.  We  will  show  that  substantial 
Improvements  may  be  obtained  by  modifying  the  boundary  conditions  of  the 
linear  stability  problem.  This  is  by  no  means  a  new  suggestion,  but  it 
seems  to  us  that  it  has  not  been  explored  in  sufficient  generality, 
especially  in  view  of  the  flexibility  which  modern  electronic  and  sensor 
technology  makes  available. 

The  plan  of  this  part  of  the  report  is  as  follows:  in  section 
1.2  we  summarize  the  relevant  aspects  of  linear  stability  theory.  In 
section  1.3  we  discuss  methods  for  manipulating  the  boundary  layer 
profile  about  which  the  linear  stability  analysis  is  done.  In  section 
1.4  we  introduce  the  general  notion  of  adaptive  boundary  conditions  and 
their  influence  on  linear  stability. 

1.2  The  Linear  Stability  Problem 

For  the  purpose  of  studying  the  drag  due  to  skin  friction  we 
may  limit  our  attention  to  essentially  parallel  flow  past  a  flat  plate. 
The  flow  field  is  then  two-dimensional  (we  choose  coordinates  x 
parallel  to  the  plate  in  the  flow  direction  and  y  normal  to  the  plate) 
and  the  relevant  variables  are  the  x  and  y  components  of  velocity 
(u  and  v)  pressure  (P)  and,  if  needed,  temperature  (T).  The  boundary 
conditions  are  usually  that  u  ■  v  ■  0  at  y-0  and  that 
u  ■  U  at  •  .  The  boundary  layer  profile  is  governed  by  a  diffusion- 
type  equation  in  which  x  plays  the  role  of  time.  This  has  the 
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consequence  Chat  Che  boundary  layer  thickness  which  we  denoCe  by 
5  ,  grows  (typically  as  /x  )  as  x  Increases.  The  actual  profile 
depends  on  the  precise  conditions  of  the  problem. 

The  linear  stability  problem  Is  posed  as  follows:  take  a  basic 
profile,  UQ(x,y),  which  varies  so  slowly  with  x  that  Its  x^dependence 
may  be  Ignored.  Then  construct  the  equation  of  motion  for  small 
perturbations,  u^(x,y),  on  that  profile  and  look  for  solutions  which  are 
periodic  In  space  and  time  (it  Is  convenient  to  use  the  stream  function 
of  the  perturbed  flow 


|4  and  vx 


3t|> 

3x 


(1.1) 


and  look  for  solutions  of  the  form  \|>(x,y)  »  $(y)elo^x“ct^  .)  The 
equation  which  emerges  is  the  Orr-Sommerfeld  equation, 


(U(y)-c)(*"-u2*)  -  U"(y)$  -  L.  (*""-2<xV'  +  a%) 


(1.2) 


where  all  velocities  are  measured  In  units  of  U  (In  particular 
U(y)  -  Uq ( y ) ]>  and  all  lengths  are  measured  in  units  of  5  and 
R  ■  U^fi/v  where  v  is  the  viscosity  of  water.  R  is  what  is  usually 
called  the  Reynolds  number  of  the  unperturbed  flow.  In  order  for  $  to 
give  an  acceptable  velocity  perturbation,  it  must  satisfy  boundary 
conditions  $  *  ■  0  at  both  y  «  0  and  y  ■  •  •  We  are  therefore 

faced  with  an  eigenvalue  problem:  for  given  R  and  a  we  must  solve 


for  e  and  its  associated  eigenfunction.  The  phase  velocity  c  will 
In  general  be  complex  and  the  background  flow  la  stable  to  saall 
perturbations  only  if  lac  <  0  for  all  a  (in  that  case,  any  saall 
disturbance  will  decay  away  exponentially).  A  typical  plot  of  atablllty 
regions  is  given  in  Fig.  1.1.  The  key  generic  feature  is  that  for 
R  <  R  »  disturbances  of  any  wavelength  decay  with  tlae  and 
therefore  the  background  flow  is  stable.  For  R  >  R  ^  ,  soae 
disturbances  grow  exponentially  with  tine  and  the  mean  flow  is  probably 
unstable  to  transition  to  turbulence.  Since  6  grows  as  we  aove 
downstream  the  criterion  for  laminar  flow  over  a  plate  of  length  L  la 
that  U^SCD/v  <  Rcrit  >  where  g(L)  is  the  boundary  layer  thickness 
at  the  downstream  end  of  the  plate.  To  enlarge  the  regime  of  laminar 
flow,  ic  is  obviously  necessary  to  make  Rcrlt  as  large  as  possible. 


*crit 


jure  1.1 
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There  are  two  methods  of  influencing  R  .  ,  both  of  which  we 

crit 

shall  discuss.  On  the  one  hand,  the  value  of  Rcrit  depends  on  tlle 
details  of  the  background  profile,  UQ(y)  ,  and  it  may  be  possible  to 
change  the  conditions  of  the  problem  so  as  to  change  U0(y)  a  waY 
which  increases  Rcr^t  •  On  the  other  hand,  one  might  take  whatever 
background  profile  comes  naturally  but  change  the  properties  of  the 
wall,  thus  changing  the  boundary  conditions  of  the  linear  stability 
eigenvalue  problem.  We  shall  see  that  it  is  possible  in  this  way  to 
greatly  change  the  stability  of  a  flow.  In  what  follows  we  will  try  to 
summarize  what  we  know  about  both  approaches. 


where  v  is  the  kinematic  viscosity  of  water  and  Z  is  the  thermal 
conductivity  per  unit  mass. 

The  above  equations  are  parabolic:  given  a  profile  at  an 
initial  value  of  x  ,  the  equations  may  be  integrated  forward  in  x  to 
produce  a  predicted  profile  at  any  downstream  x  .  Consequently,  to 
specify  the  boundary  layer  one  must  specify  an  initial  profile — i.e. , 
one  must  know  the  details  of  how  it  started.  This  annoyance  may  be 
circumvented  by  looking  for  self-similar  flows — flows  in  which  the 
profiles  at  different  downstream  stations  differ  only  in  the  value  of  a 
transverse  dimensional  scale.  The  Reynolds  scaling  property  of  these 
equations  guarantee  that  we  can  find  solutions  in  which  the  stream 
function,  f  ,  and  temperature,  T  ,  have  the  form 


i|i(x,y)  -  /vxU  f(n) 

00 

T(x,y)  -  t(n) 

fa 

n  ■  ■  (!•*) 

The  partial  differential  equations  for  the  boundary  layer  then  collapse 
to  non-linear  ordinary  differential  equations  whose  solution  is 
completely  determined  by  the  boundary  conditions  at  y  -  0  and  y  ■  <•  . 
We  will  shortly  discuss  the  properties  of  these  solutions.  There 
remains  the  question  why  we  are  Interested  in  self-similar  solutions: 

How  do  we  know  that  the  initial  profile  is  Just  right  to  set  up  the 
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self-similar  boundary  layer?  We  will  eventually  show  that  the  self¬ 
similar  solutions  have  a  certain  stability  property  which  guarantees 
that  small  deviations  from  the  self-similar  profile  die  away  as  we  go 
downstream.  We  take  this  as  a  Justification  for  concentrating  on  self¬ 
similar  profiles  for  our  discussion  of  the  linear  stability  problem 
outlined  in  the  previous  section. 

Let  us  first  summarize  the  properties  of  the  standard  flat 
plate  self-similar  boundary  layer  in  which  temperature  variations  are 
ignored  and  the  standard  wall  boundary  conditions  (u  =  v  =  0  at  y  ■  0) 
are  applied.  The  equations  to  be  solved  are 


2f "  '  +  ff"  «  0  (1.5) 

i 

u  =  u^f'Cn)  ,  v  -  (nf'-f)  (1.6) 

f=f'=0  at  n  0  (1.7) 

f '  =  1  at  n  =*  °°  (1.8) 


The  equation  for  f  must  be  solved  numerically. 

We  need  a  definition  of  boundary  layer  thickness  to  specify  a 
Reynolds  number.  There  are  two  physically  significant  definitions: 
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Displacement  Thickness  -  -  f  dy(i-  ^-) 


(1.9) 


Momentum  Thickness  •  £_ 


■  /  S-  (i-  5-) 

J  mm 


(1.10) 


iSj  is  the  net  outward  displacement  of  a  far  field  streamline  due  to 
the  presence  of  the  plate.  6^  measures  the  momentum  lost  from  the 
flow  due  to  the  presence  of  the  plate.  For  the  self-similar  flow  at 
hand 


/TP  [n-f(n)J  rr(B 


1.72 


f  \K 

IV~ 


(1.11) 


6 


2 


(1.12) 


Both  thicknesses  grow  as  /x  as  we  move  downstream  from  the  front  edge 
of  the  plate.  It  is  customary  to  use  5^  to  construct  the  Reynolds 
number  for  the  flow: 


U.6i 


(1.13) 


The  ratio  6 ^ »  which  in  this  case  has  the  value  2.59,  turns  out  to 
be  a  convenient  diagnostic  quantity  for  stability  of  the  flow. 


The  classic  linear  stability  theory,  when  applied  to  this  flow, 
gives  a  critical  Reynolds  number 


crit 


420 


(1.14) 
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This  Is  rather  small — If  we  want  absolutely  stable  flow  over  a  surface  1 

m  long,  the  maximum  flow  velocity  Is  6  cm/s!  This  Is  due  to  the  fact 

that  Blaslus  flow  is  just  on  the  boundary  of  instability.  If  the 

2  2 

velocity  profile  has  a  point  of  Inflection  3  U/3y  -  0  ,  the  flow  is 
potentially  unstable;  Blaslus  flow  has  a  point  of  Inflection  precisely 
at  y  -  0.  To  increase  stability  it  is  obviously  necessary  to  decrease 
the  curvature  at  y  *  0  by  some  means. 


are  now  constants 


(1.16) 


and  R  ■  U  /v  anywhere  along  the  plate.  The  linear  stability 
Oj  •  o 

analysis  shows  that  Rcr|t  “  50,000  ,  a  factor  10^  improvement  over  the 
no  suction  value.  This  is  in  part  due  to  the  fact  that  the  curvature  of 
the  velocity  profile  is  now  everywhere  negative  and  in  part  due  to  the 
fact  that  the  diagnostic  quantity  has  been  reduced  to  2  from 

1-9 
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2.6.  The  practical  possibilities  are  Impressive:  If  one  chooses 
U  /v  -  10^  ,  well  within  the  stability  region,  It  Is  possible  to 

oe  o 

stabilize  the  boundary  layer  of  a  5  m  diameter  100  m  long  submarine 
moving  at  20  kts  by  ingesting  water  at  the  trivial  net  rate  of  3  m^/s! 
The  problem  Is  that  In  order  for  this  analysis  to  apply  the  suction  must 
be  applied  uniformly  over  transverse  dimensions  smaller  than  the 
boundary  layer  thickness  Itself.  Since  the  boundary  layer  Is  only  a 
fraction  of  a  millimeter  thick  this  means  that  the  boundary  layer  has  to 
be  Ingested  through  small  pores  which  will  almost  certainly  clog  in  real 
ocean  use.  An  alternate  approach  which  may  overcome  this  problem  will 
be  discussed  at  the  end  of  this  section. 

The  second  stabilization  method  is  to  Increase  the  temperature 
of  the  wall  above  that  of  the  mean  flow.  Since  the  viscosity  of  water 
at  room  temperature  decreases  rather  rapidly  with  Increasing 
temperature,  It  becomes  Important  to  include  the  temperature  in  the 
boundary  layer  equations.  If  we  evaluate  the  first  boundary  layer 
equation  at  the  wall  (where  u  ■  v  *  0) ,  we  find: 


If  the  wall  Is  hot,  viscosity  must  Increase  outward  as  must  u  . 

9  9 

Therefore  (3*U/3y*)g  Is  less  than  zero,  a  condition  we  have  already 
argued  must  Increase  stability. 
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Numerical  Integration  of  these  equations  for  a  heated  wall  have 

been  carried  out  by  several  authors.  It  is  found  that  (R,  1  ... 

v  'crit 

Increases  dramatically  with  Increasing  heating,  with  a  maximum  value  of 

I 

12,000  being  obtained  for  50°C  temperature  difference.  This  great 
Increase  In  stability  Is  not  achieved  at  the  cost  of  extreme 
modifications  of  the  profile.  What  seems  to  count  Is  the  ratio 
6^/62  '•  Rcrlt  correlates  very  well  with  this  quantity,  Increasing  by 
10^  as  decreases  from  2.6  to  2! 


Finally,  we  want  to  return  to  the  problem  of  using  suction  to 
control  the  boundary  layer.  An  alternative  to  distributed  suction  Is  to 
extract  the  boundary  layer  fluid  through  large,  widely-spaced  slots. 

The  idea  is  that  the  boundary  layer  Is  allowed  to  grow  for  a  distance 
D  and  then  is  sucked  Into  an  internal  reservoir  through  a  slot  large 
enough  to  avoid  clogging.  A  fresh  boundary  layer  forms  at  the 
downstream  edge  of  the  slot  and,  after  a  distance  D  ,  Is  sucked  in 
again.  The  maximum  thickness  of  the  boundary  layer  Is  evidently 


6 


max 


(1.18) 


In  order  for  that  thickness  of  fluid  to  be  sucked  In  as  it  passes  a  slot 
of  width  d  at  stream  velocity  U  ,  we  must  apply  a  suction  velocity 


v  7  jv_  d  1  d 

U  n/UD  D  ~  R_  D 

5 

max 


(1.19) 


The  condition  for  stable  flow  along  the  plates  of  length  D  Is  that 
Rg  <  400  .  This  Is  achieved  for  Interesting  values  of  for 

D  -  1  m.  We  could  realistically  choose  d  -  1  cm  and  therefore  control 
the  boundary  layer  with  ^  <  10”^  ,  which  corresponds  to  a  very 
manageable  rate  of  water  Ingestion*  The  difficulty  might  be  Instability 
associated  with  the  new  type  of  shear  layer  across  the  slot  Itself*  It 
Is  not  clear  that  this  shear  layer  Is  any  more  unstable  than  that  along 
the  wall,  but  we  have  seen  no  discussion  of  this  problem.  To  repeat, 
the  essence  of  this  method  is  to  achieve  stability  by  interrupting  the 
growth  of  the  boundary  layer  every  D  meters  so  that  the  effective 
scale  size  for  the  stability  problem  Is  D  rather  than  the  overall 
vehicle  length* 

1.4  Adaptive  Boundary  Conditions 

The  discussion  so  far  has  dealt  with  situations  In  which  the 
boundary  conditions  at  the  wall  are  fixed,  and  we  have  seen  that  large 
modifications  In  stability  can  be  obtained  with  "small"  changes  in  those 
boundary  conditions*  It  might  be  more  efficient  to  modify  the  wall 
conditions  only  in  response  to  disturbances  on  the  basic  flow,  trying  to 
tailor  the  response  In  such  a  way  as  to  damp  the  disturbance.  If  we  are 
very  lucky  a  scheme  of  this  kind  might  allow  us  to  obtain  higher  values 
of  Rcrlt  * 


The  basic  situation  is  as  described  in  section  1.2:  We  have  a 
background  flow  described  by  a  horizontal  velocity  profile  UQ(y)  . 
Small  perturbations  on  this  flow  are  described  by  a  stream  function 
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i|i(x,y)  »  $(y )eia^x~ct)  which  satisfies  the  Orr-Sommerfeld  equation. 

We  may  eventually  want  to  expand  this  system  to  Include  the  effect  of 
temperature  variations  or  to  allow  the  background  profile  to  have  inflow 
due  to  suction;  but  we  will  stick  with  the  simplest  system  for 
explanatory  purposes. 


J 


We  propose  to  modify  the  nature  of  the  eigenvalue  problem  for 
c  by  Imposing  a  non-standard  boundary  condition  at  the  wall  (we  of 
course  retain  the  usual  zero  velocity  condition,  $  »  »  0,  at  y  *  »  ) 

In  which  the  flow  at  the  wall  Is  a  linear  function  of  the  past  history 
of  flow  quantities  sensed  at  the  wall.  Let  S(x,t)  be  the  quantity 
sensed  at  the  wail  and  let  y(x,t)  be  the  response.  Then  the  sort  of 
boundary  condition  we  are  contemplating  is 


Y(x,t) 


dyG(x-y,t-t' )S(y  ,t 


)  . 


(1.20) 


Since  the  small  disturbances  are  taken  to  depend  on  x  and  t  as 
ela(x-ct)  ^  the  above  equation  may  be  Fourier  transformed: 


y(a,ac)  -  ?(a,ac)§'(a  ,ac)  •  (1.21) 

The  only  essential  property  of  d  Is  that  the  singularities  in  ac  lie 
In  the  lower  half  plane  so  that  G  is  causal.  G  may  reflect  the 
mechanical  properties  of  a  material  surface,  in  which  case  C  Is 
prescribed,  or  it  may  be  constructed  via  electronic  coupling  of  sources 
and  drivers,  in  which  case  it  may  be  whatever  we  like.  The  general 
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"  -J 


problem  will  be  to  use  the  freedom  to  pick  S  and  G  to  optimize 
stability.  Subsets  of  this  problem  have  been  analyzed  In  some  detail, 
but  an  overview  seems  to  be  lacking. 


Let  us  summarize  some  of  our  options  for  sensing  and 
responding.  In  terms  of  the  perturbed  stream  function, 
m  $(y)ela(x-ct)  ,  the  field  quantities  are 


u 


v 


lap 

~P~ 


*'(y)ela<x-ct> 

(1.22) 

-la*(y)ela<x-ct> 

(1.23) 

v  /-a^+  — +  iaU'$  - 

ia(U-c)<|>' 

ia(x-ct) 

e 

\  dy  / 

- 

(1.24) 

The  wall  Is  subject  to  transverse  forces  (due  to  pressure)  or 
longitudinal  forces  (due  to  Reynolds  stress,  t  »  p(-|^  +  -|^)  )  and 

either  could  be  the  quantity  we  want  to  sense.  The  sensed  quantity  will 

therefore  be  some  linear  combination  of  $  and  its  derivatives 

evaluated  at  the  wall.  The  response  quantity  could  be  a  motion  of  the 

wall  Itself  (compliant  wall)  in  either  the  normal  or  the  transverse 
direction  or  a  suction  velocity  without  physical  wall  motion  (adaptive 
suction).  Again,  the  response  quantity  will  be  a  linear  combination  of 
t  and  its  derivatives  evaluated  at  the  wall.  If  the  wall  moves  as 
part  of  its  response,  the  wall  is  no  longer  at  y  ■  0,  a  fact  which  must 
be  accounted  for  in  writing  out  the  boundary  condition  on  $  .  Clearly, 
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temperature  or  heat  flux  could  be  added  in  an  appropriate  way  to  this 


list. 


We  now  list  some  specific  boundary  conditions  of  the  type  we 
are  discussing.  In  all  cases  the  problem  is  to  solve  the  0-S  equation 
in  the  interval  0  <  y  <  «  with  $(«)  ”  $’(“>)  =  0  plus  the  two 
specified  conditions  at  y  *  0: 

Normal  Compliant  f*,vll: 

c<J> '  ( 0)  +  U'(0)$(0)  =  0  ;  io*(0)  =  G(a.ac) 

(1.25) 


Transverse  Compliant  Wall: 

4>'(0)  =■  C'(a,ac)4."(0)  ;  4,(0)  =  0  (1.26) 

Adaptive  Suction  (sensing  pressure) 

(0)  *  0  ;  ia$(0)  *>  (T(a,ac)  ^  —  (1.27) 

P 

Adaptive  Suction  (sensing  Reynolds  stress): 

$'(0)  -  0  ;  ia*(0)  *  G'(a,ac)$' '  (0)  (1.28) 
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In  all  these  cases  the  first  condition  refers  to  x  components  of 
velocity  and  the  second  to  y  components*  p(0)  is  meant  to  be 
evaluated  in  terms  of  $  .  As  advertised*  the  new  boundary  conditions 
Involve  linear  combination  of  $  and  its  derivatives,  plus  an 
adjustable  response  Green's  function  which  we  have  denoted  by  ST  . 
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THE  EFFECT  OF  BOUNDARY  CONDITIONS  ON  STABILITY  OF  FLOWS 


2.1  Introduction 

Many  suggestions  have  been  made  to  reduce  the  drag  on  bodies 
moving  through  fluids  by  delaying  the  transition  from  laminar  to 
turbulent  flow.  Generlcally  these  Ideas  can  be  divided  into  two 
classes.  One  Is  to  modify  the  flow  In  the  vicinity  of  the  body. 

Examples  are  heating,  sucking,  and  the  dispersal  of  polymers.  The  other 
is  to  change  the  nature  of  the  surface  presented  by  the  body  to  the 
fluid.  This  can  be  done  either  actively  or  passively.  Here  we  ask  what 
can  be  said,  very  generally  about  the  passive  case — compliant  boundary 
conditions. 

As  one  of  the  simplest  possible  models  we  consider  the  case  of 
invlscld,  Incompressible  parallel  shear  flow.  We  imagine  a  flow  U(z) 
established  between  plates  at  Zj  and  z2  (z2  may  be  infinity).  The 
stability  of  the  flow  Is  then  examined  as  a  function  of  the  boundary 
properties. 

2.2  The  Classical  Theorems 

The  case  In  which  the  boundaries  are  both  rigid  plates  has  been 
extensively  studied. 

By  virtue  of  Squire's  theorem  we  can  restrict  attention  to 
purely  two  dimensional  perturbations.  Thus  we  consider  a  perturbed 
stream  function  of  the  form 


II-l 


*  *  $(z)e 


i  (  OK-C  t  ) 


(2.1) 


The  stability  problem  reduces  to  finding  (for  fixed  real  a  )  when  there 
are  discrete  eigenvalues  c  ■  cr  +  ic^  with  c^  >  0  of  the  Rayleigh 
stability  equation 

[*l  '  *2)*  *  *  0  (2-2> 


with  the  boundary  conditions 


*(Zl)  -  *(z2)  -  0  (2.3) 

Fortunately  there  are  some  (not  many)  general  results  known 
about  this  problem.  These  include: 

(1)  Rayleigh's  Inflection  Point  Theorem. 

This  is:  A  necessary  condition  for  Instability  is  that 
the  basic  velocity  profile  should  have  an  inflection  point,  i.e.,  there 
is  a 


z  ,z  «  z  <  z  such  that  U''fz  )  ■  0  (2.4) 

s  1  s  2  v  s 

(2)  Fjortofta'  Theorem 

A  necessary  condition  for  instability  is  that 


U"(U-U  )  <  0  (2.5) 

somewhere  in  the  flow.  (Here  U^s  U(z#)  ]. 
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( 3 )  Howard's  Semi -Circle  Theorem 


For  unstable  waves  c  must  lie  in  the  semi-circle  with 

radius  ifu  -  U  .  1  ,  center  at  i-fu  +  U  .  1  ,  and  c.  >  0  ,  i.e. , 

21-  max  rain-'  *  max  rain 1  *  i  ’ 

2 

r  1  \  i2  ,  2  ,  1  |  max  mini  ..  n 

c  --x-lU  +uJjl  t - 1  ,  c.  >  0 

1  r  2V  max  rain-'-1  i  2  \  2  /  i 


(2.6) 


It  should  be  emphasized  that  the  Rayleigh  Inflection  Point 
Theorem  is  a  necessary  but  not  sufficient  condition  for  instability.  It 
is,  however,  of  great  qualitative  importance.  Even  when  considering 
viscosity  most  attempts  to  delay  (or  elimina'e)  instability  are  devoted 
to  getting  as  far  away  as  possible  from  inflection  points. 


There  are,  however,  a  class  of  flows  (which  include  most 
boundary  layer  flows)  for  which  there  also  exists  a  sufficient  condition 
for  instability.  Thus  if 


K(z)  =  U"/(U-UgJ 
is  regular  at  z  =  zg  and 

K(z)  >  ■n2/(z2~z1  )2 


everywhere,  then  we  have  instability. 


(2.7) 


(2.8) 
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2.3 


More  General  Boundaries 


In  Che  case  of  compliant  boundaries  the  boundary  conditions  of 
Eq.  (2.3)  will  no  longer  hold.  The  stability  Eq.  (2.2)  will  still  be 
valid  and  since  it  is  a  second  order  differential  equation  we  can  put 
conditions  at  Zj  and  '  As  a  rather  general  case  we  take  these 
conditions  to  be  that  a  ’kiiear  combination  of  $  and  its  normal 
derivative  be  zero. 

Thus 


-  +  'l. 


(2.9) 


(Note:  The  classical  case  is 


-►<»). 


0. 


A  free  boundary  corresponds  to 


Now  let  us  see  what  happens  to  the  previous  theorems. 


( 1 )  The  Inflection  Point  Theorem 

Multiplying  Eq.  (2.2)  by  ,  integrating  from  Zj  ,  to 
Z2  and  using  the  boundary  conditions,  Eq.  (2.9)  gives 


(2.10) 
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Let 


*1  ■  si + "i  • 


Then  the  Imaginary  part  of  Eq.(2.10)  becomes 


*z2 


1 3*l 

+  T  1^1 

fazl  ! 

2\Hz\ 

(2.11) 


We  note  that  in  gene ral  the  requirement 


/' 


Cj  >  o  *  t  ~2  u- V*z  ■  o 


|U-c| 


'1 


only  holds  if  Ti  =  T2  =  0  .  That  is  the  Inflection  Point  Theorem  holds 
only  if  Y^  i  is  real.  In  this  case 


c  =  0  ♦  U"(z  )  =  0  . 

i  s 


Until  it  is  indicated  differently  we  will  be  assuming  Y^  is 
indeed  real  so  that  the  Inflection  Point  Theorem  holds. 


(2)  Fjortofts'  Theorem 

The  real  part  of  Eq.  (2.10)  is 

f2  i  a  A 1 2  ,  2  U"(U-c  )H>| 


dz 


s  |3*|2  +  s  liil2 

S1  \lt\l  +  ^ 2 1  D s'- 1  2 


(2.12) 
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Adding  to  this 


0 


(2.13) 


gives : 


U"(U-Ug)  |«|‘ 


I U— c 1 ' 


dz 


(2. 14) 


Thus  in  general  Fjortoft's  Theorem  holds  only  if  Sj  and  S2  <  0  •  Then 
a  necessary  condition  for  instability  is 

U"(U-U  )  <  0  somewhere  in  the  flow, 
s 

(3)  Howard's  Circle  Theorem 

This  theorem  essentially  completely  disappears  when  our 
more  general  boundary  conditions  are  applied.  Even  the  weaker  Rayleigh 
theorem  that  cr  must  be  in  the  range  of  U(z)  fails  to  hold.  To  see 
this  we  parallel  the  usual  derivation  of  the  theorem. 


Let 


fc+  - 


(2.15) 
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Then  Eq.  (2*12)  becomes 


O^)2  4-  /  “  a2(U-c)2  $+  -  0 


(2.16) 


Multiply  this  equation  by  $*+  and  integrating  yields 


(U-c)  dz 


rz2 

l  + 

2  2 

/ 

I  3z 

+  a2U+l 

Z  1 

[<“-)2**+^rc 


(2.17) 


Now  If  ^(Zj)  -  <p(z 2 )  -  0  then  ^(z^)  -  <f>+(z2)  =  0.  (We  are  assuming 
Ci  >  0  and  hence  U-c  j4  0  • )  The  imaginary  part  of  Eq.  (2.17)  becomes 


2'i/2|l#r*-vr 


(U-cr )dz  *  0  ,  (2.18) 


1 


and  hence  cr  is  in  the  range  of  U  .  Consider,  however,  the  general 
boundary  condition 


I 


I 


«(z1)  -  St 


i 


1,2 


Then  (2.17)  becomes 


3$ 

JT~ 


a  VI 


(U-c)  dz 


(U-c) 


1  U' ' 
'  tPc 


(U-c; 
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1  +  U' 
5^  +  U=c 


UT 


(2.19) 


(2.20) 


Assuming  ci  t  0  we  obtain  from  the  imaginary  part  of  this  equation  that 


P  j|£|2  +  aVlVt,d,- 

z. 


■ 

L  >  ifi2l 

4. 

(U-cr)U+|2 

Z2 

S1 

*1 

Iir 


u+r 


(2.21) 


No  particularly  interesting  result  seems  to  follow  from  this. 


2 •  A  Positive  Results 

So  far  all  we  have  seen  is  that  the  introduction  of  our  general 
boundary  conditions  weakens  the  few  general  theorems  known  about 
parallel  shear  flows.  Now  we  will  show  that  paralleling  the  proof  of 
the  sufficient  condition  for  instability  gives  us  a  method  to  estimate 
the  stability  dependence  on  the 

Tbe  essence  of  the  approach  is  the  following:  We  construct  a 

neutrally  stable  solution  for  a  particular  a2.  Then  one  shows  that 

8 

for  a  slightly  lower  than  o2  there  is  an  unstable  solution.  Let 
K(z)  -  -U"/(U-U8)  be  regular  at  zg  .  Set  c  •  c8  •  U8» 

The  Rayleigh  Equation  (2.2)  becomes 

♦"  +  [K<*)-a2  ]♦  -  0  (2.22) 
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This  Is  a  Sturm-Liouville  problem  with  the  variational  principle 


/V I2-Kum,v  - 


First  we  show  how  this  can  be  used  to  obtain  a  sufficient  condition  for 
instability  and  then  how  it  determines  dependence  on  the  S, „ 


Varying 

f  we  see  that  the  minimum  is 

obtained  if 

f "  +  nf  *  0 

and 

fUj)  -  -  st  f'(2l) 

f(z2)  -  S2f'(z2) 

Let  n 

be  the  smallest  eigenvalue  of 

this  problem 

therefore, 


< 


K(z)|$|  dz 


M>|2dz 


(2.26) 


(2.27) 


Thus,  if  K(z)  >  n  everywhere 
m 


2 


a 


s 


<  0 


and  so  there  is  a  neutrally  stable  mode.  (Below  we  will  calculate 
for  the  special  case  S2  “  0  •  If  will  be  shown  that 


^m 


(2.28) 


where 


0  <  f(Sj )  <  1  . 


(2.29) 
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t 


I 


I 


I 


It  can  be  shown  that  if  there  is  such  a  neutrally  stable  mode, 

2  2 

then  there  is  an  unstable  mode  for  a  Just  less  than  ag  •  This  can 
be  proved  either  by  using  perturbation  theory^ ^  or  more  conveniently 
using  the  variational  principle  as  demonstrated  in  the  Appendix. 


We  can  now  see  the  effects  of  varying  the  on  stability. 

Suppose  we  have  for  some  values  of  Sj  ,  S2  a  neutrally  stable  mode 

(and  consequently  an  unstable  mode  nearby).  Then  if  by  changing  the 
2 

we  can  force  ag  to  become  negative  we  will  achieve  stability.  But 

2 

from  Eq.  (2.23)  we  readily  find  the  change  of  ag  with  the  . 

2  2 
Indeed  since  -  ag  is  stationary  we  obtain  the  rate  of  change  of  ag 

merely  by  differentiation  with  respect  to  the  explicit  dependence  on 

the  .  Thus 


(2.3u) 


and 


f 


|$(z)|  dz 


(2.31) 


We  see  these  slopes  are  always  positive.  Thus  any  Increase  in  Sj  or 

2 

S2  always  pushes  -ag  to  zero.  The  most  stable  situation  is  then  that 
with  Sj,  Sj  ♦  »  .  The  boundary  conditions  are  then  "  0,  i.e.,  a 

free  boundary. 


II-l  i 


*•  V 


There  is  a  caveat  here,  however.  In  general,  the  dispersion 

a 

relation  lor(c)]  has  many  branches.  Since  we  have  implicitly  assimed 
a  continuous  behavior,  what  has  really  been  proved  is  that  corresponding 

2 

to  a  given  branch  -  ag  always  increased  with  increasing  Sj,S2  •  It 
may  happen  that  corresponding  to  some  values  of  the  St  new  branches  of 
the  dispersion  relation  may  appear — and  a  new  instability  may  arise. 

All  this  is  exemplified  by  the  following  example. 


A  Classic  Example 

Consider  U  -  sin  z  ,  z.  <  0  <  z„ 

1  2 


Then  there  is  an  inflection  point  at  zg  -  0  ,  and  „e  can  take 
cs  “  0  .  The  stability  equation  becomes 


♦"  +  (!-«:?)♦  -  o 


(2.32) 


For  simplicity  we  will  choose  S2  -  Q  ,  i.e.,  «,(z,)  -  0  . 


(2.33) 


(Note:  The  n  of  the  previous  section  is  Just  1  -  ag  .  )  Using  the 
boundary  condition  Eq.  (2.33)  and 


♦(*i)  ■  -s 


(2.34) 


we  find 


“a2"1 


(2.35) 
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I 


where  X  is  a  solution  of  the  equation 

tan  X  -  - - -  X  •  (2.36) 

z2  ~  Z1 

Thus  X  is  the  intersection  of  the  curve  y  »  tan  X  with  a  straight 
line  through  the  origin  with  slope  S^/Cz^-  z^). 

S1 

Let  us  start  with - -  <*>  . 

z2~zi 

The  intersections  occur  at  (n  +  1/2 )  tt,  n  »  0  ,  1  ,  .  .  . 

F  rom  Eq .  (2.35) 


(n+l/2)ir]2 

Zl  "  Z2j 

and  we  have  instability  only  if 

(n+l/2)u  <  (z2_z1) 


(2.37) 


(2.38) 


Clearly  the  mode  n  »  0  is  the  most  unstable  and  the  requirement  for 
stability  is 


z2  "  zt  <  */2 


As  Sj/(z2~Zj)  increases  from  •  «  to  0  the  points  of  intersection  move 
smoothly  from  (n  +  1  /  2 )  it  to  (n  +  1 )  it  .  The  n  -  0  mode  is  still  the 
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most  unstable.  At  Sj  ■  0  we  have  stability  If  Zj  “  <  »  •  Thus 

the  range  of  Z2  -  zj  for  which  we  have  stability  has  doubled. 

As  sj(z2  ”  zj)  grows  from  0  to  +  »  intersection  points  move 

from  (n  +  l)ir  to  (n  +  3/2 ) w  .  Again  the  most  unstable  of  these  Is 

3 

n  »•  0  and  we  apparently  have  stability  for  z^  -  z^  <  .  •  This  Is, 

however,  Incorrect.  For  S^/fz^-  z^)  ■  1  a  new  mode  appears.  Thus  with 
X  ■  0  we  have  a  solution 

$  *  z^  -  z  .  (2.39) 

This  Is  quite  unstable  with 

a2  “  1  .  (2.40) 

s 

When  S^/Zj-  zp  Increases  from  l  to  +  •  the  Intersection  for 
this  new  mode  moves  from  0  to  ir/2  .  Thus  In  this  range  we  have 
stability  only  for  (zj  -  z^)  <  and  for  part  of  the  region 

z2  ~  Z1  <<  2  * 

Combining  the  results  we  see  that  the  most  stability  Is 
obtained  when  S^Zj  -  z^)  ■  1  -  e  .  Then  we  have  stability  for 

(z 2  -  2j)  <  w(l  +  6)  (2.41) 
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where 


0  <  5  < 


6-1/4  . 


(2.42) 


2.6  Complex  Boundary  Conditions 

In  the  above  we  have  found  that  for  real  boundary  conditions 
some  general  global  statements  can  be  made*  Thus,  the  Inflection 
Theorem  holds,  Fjortoft's  theorem  holds  for  <  0,  ,  and  for  a  given 
branch  of  the  dispersion  curve  increasing  the  increases  stability. 


With  complex  boundary  conditions  we  have  not  found  any  global 
theorems.  However,  we  can  make  a  local  statement.  Thus  in  the  Appendix 
it  is  shown  that  in  the  vicinity  of  a  neutrally  stable  solution  there 
are  two  complex  conjugate  solutions  for  o  <  as»  One  of  these  is 

unstable.  By  an  appropriate  choice  of  complex  it  is  possible  to 

delay  the  instability. 


We  have  the  expressions 


3c  _  ,  . 

15  P1  i  1  p2 


(2.43) 


C2 1,^.2  Au"  2i j  i*m2 

J  !»♦  *  +  1*1  ldz  “  - V - - v-1- 


2  *1  ' 


l$(z) I  dz 


(2.44) 


» 
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The  change  in  -a2  at  -a2  when  y,  -  S.  is  replaced  by  s.  +  16T, 

ail  l  i 

is  (since  a2  is  stationary)  obtained  merely  from  the  explicit  change 
s 

in  the  Yj  l*e* » 


6(-a2)  -  iy  i5T2p(z2)|2 
S2 


+  1  1ST  |*(z  )|2 

S1 

i.e.,  6a  “  iB  with  B  real. 

s 


(2.45) 


Then 


6c  -  (pj  +  ip2)iB  -  iBp^  +  Bp2  . 

With  the  appropriate  choice  of  sign  of  B  we  find 

Im  6c  <  0. 

Therefore,  locally  a  complex  Y^  (phase  lag)  can  be  chosen  to 
stabilize. 


(2.46) 


(2.47) 


2.7  Conclusion 

It  has  been  found  that  for  real  boundary  conditions  the 
stability  for  a  given  branch  of  the  dispersion  relation  is  always 
Increased  by  increasing  S|  (i.e.,  the  closer  to  a  free  boundary  one  is 
the  greater  the  stability).  There  is  a  warning,  however.  New,  more 
unstable  branches  can  arise. 
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t 


Using  complex  boundary  conditions  can  at  least  locally  delay 
>  Instability. 

Hopefully  in  a  following  note  we  will  endeavor  to  see  what  can 
)  be  said  when  viscosity  is  included. 

< 


t 
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COMPLEX  BOUNDARY  CONDITIONS 


1  Introduction 

We  have  seen  that  when  the  classical  rigid  body  boundary 
conditions  are  replaced  by  more  general  ones  very  little  can  be  said  in 
generality.  Accordingly  it  is  useful  to  consider  some  simple  models. 

This  is  done  here.  It  is  found  that  dramatic  changes  in  stability  can 
be  achieved. 


The  situation  we  envisage  is  that,  instead  of  vanishing  at  the 
boundary  (z  -  Zi),  the  stream  function  should  satisfy 

£  (x-.ZjMx-  .  <3.1; 

—00 

When  Fourier  transformed  in  the  x  coordinate  this  becomes 


► 


1 


> 


*(zi>  *  v*)  u  <2i)  • 

The  reality  condition  on  G  implies 

V~a)  *  • 

Alternatively  this  can  be  written: 


Ill-i 


(3.2) 


(3.3) 


Let 


Y1(a)  -  S^a)  +  iT^(a) 


(3.4) 


The  reality  condition  is  then  that 

Sl(-a)  -  St(a) 


(3.5) 


and 


Tt(-a)  -  -  Tt(a).  (3.6) 

Since  the  Rayleigh  stability  equation  is  rather  intractable  for 
a  general  parallel  shear  flow  we  restrict  ourselves  to  models  in  which 
the  basic  flow  is  at  most  piece-wise  linear  in  z  .  However, 
discontinuities  in  U  and  its  derivatives  are  admitted. 


The  problem  then  is  to  find  those  c  such  that 


a  $ 


subject  to  the  conditions 


(3.7) 


M(U-c)$'  -  U'*]  -  0  , 


(3.8) 


4 


Mel  -  0  • 


(3.9) 


at  discontinuities  and 


1 1 1-2 


» 


♦  C*i)  -  #  C*i  )  .  *-1.2  »•“» 


at  boundaries. 


3.2  Model  I 


We  consider  the  flow 


U  =  0  ,  0  <  z  <  d 


(3.11) 


U  =  Sz 


z  >  d 


(3.12) 


and  the  boundary  conditions 


♦(0)  -  -Y(a)  ^  (0),  ♦(«)  »  0  .  (3.13) 


The  eigenvalue  condition  then  is 


clal  =  ~  1  - 


.2  (1  +  |o|Y)| 
e  (1  -  lafTT) 


(3.14) 


where  e  =*  exp(-ad). 


In  particular  for  the  imaginary  part  of  c  we  obtain 


c^  »  -  Se 


2  Im  Y 
[  1-1  a|Y 


(3.15) 


To  incorporate  the  reality  properties  we  write 


1 1 1-3 


I 


(3.16) 


Y  -i^“>  +  i  r<°> 
I  a  I  o 


Then  S'  and  T'  must  be  even  functions  of  a  •  The  real  growth  rate 
ac^  Is  then  given  as 


-  Se2T' 

ac .  “  - 5 - 7— 

[a-s')  +  r  } 

If  we  maximize  with  respect  to  T'  we  find 


(3.17) 


T'  -  1  -  S' 


(3.18) 


and  then 


ttCl  -  •  (3. 19) 

Letting  S'  go  to  1  from  below  or  above  (depending  on  the  sign  of  S  ) 
we  see  that  for  all  a  we  can  obtain  arbitrarily  large  damping.  The 
corresponding  Y  Is 


Y  -  -  I  +  1  9&nS  n  (3.20) 

I  a  I  a 

with  n  arbitrarily  small  and  positive. 

3.3  Model  2 

This  Is  similar  to  the  above  but  slightly  more  complicated.  We 
take  boundaries  at  Z}  and  Z2  with  the  same  Y  , 
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♦<*1>  ■  Y<°>  u  (*«> 


i  -  1,  2.  (3.21) 


The  unperturbed  flow  is 


U  -  Uj, 


Zj  <  z  <  d  +  z  ^  (3.22) 


U  -  U  +  S(z-z  j-d), 


z^+  d  <  z  <  z2  -  d  (3.23) 


U  =  U2  “  U1  +  S(Z2~  d)’  z  >  Z2  "  d  *  (3-24) 


Since  the  flow  is  so  symmetric  we  can  divide  the  solutions  into  even  and 
odd  ones  with  respect  to  the  midpoint  which  we  take  to  be  z  =■  0  .  Thus 
we  merely  need  to  look  for  solutions  for  o  <  z  <  subject  to 
<f>(0)  =  0  or  <J>'(0)  *  0  .  Also  without  loss  of  generality  we  can  take 

U2  =  0.  We  choose  units  so  z2  =  1  . 

The  eigenvalue  equation  then  becomes 


TToT 


J  jl-  -i2-j[l+la|Y]e-2|a|-[e2-x](l-|a|Y]j 

- -  _ -  (3. 

n  j  i  ..  ivu  4  i  \  \  _  ,  %r  \ 


25) 


[ 1+1 o|Y]e  -  X(l-|a|Y) 


Here  c  is  again  exp  -  |a|d  and  \  *  +1  for  odd  solutions  and  -1 
for  even  ones. 

This  is  somewhat  complicated  to  discuss.  Let  us,  however,  just 
look  at  the  region  la|  »  1,  d  ~  1  . 
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The  equation  then  becomes 


C  - 

,  4.  »-2a(l-d)  [  1  +  1  a| Y J 

7a 

1 1-1  at  YJ  J 

(3.26) 


This  then  is  just  as  for  our  first  model.  Therefore  again  we  can 
achieve  arbitrarily  large  damping. 


3* 4  Model  3:  The  Helmholtz  Instability 

Consider  the  flow  in  the  region  z^  <  z  <  z2  where  Zj  <  0  , 
z2  >  0  .  For  Zj  <  z  <  0  ,  U(z)  -  -UQ/2  and  for  z2  >  z  >  0  , 

U(z)  -  Uq/2.  As  boundary  conditions  we  take 

♦  (*!)  -  0  •  *(*2)  “  Y  ~ g~  •  <3*27> 

If  x  -  [Uq/2-c ]/[Uq/2+c ]  the  eigenvalue  equation  becomes 

2  tanhf | ol z2  ]  -  |a|Y 

*  - - - - — -  (3.28) 

jtanh[ laiZj  ]  }{l-|a|Y  tanh|a|z2) 

Note:  If  x  is  found,  the  c  Is  given  by 


c 


u0  (1  -  x) 
7"  (I  +  x)  * 


(3.29) 


First  consider  Y  =  0  (the  usual  Helmholtz  case.)  Then 

2  tanh[|a|z2] 

X  •  1  ■■  • 

tanh[ |o|«1 ] 


(3.30) 
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Since  Z2  >  0  ,  and  zj  <  0  we  see  x^  is  negative  for  all  a  • 
two  roots  for  x  are  purely  imaginary  and  so  one  of  them  corresponds  to 
an  unstable  mode. 

Suppose  now  Y  is  real  and 

I a| Y  ♦  1 /tanh[ | a| z^ ]  .  (3.31) 

Since  tanh^(|a|Z2)  <1  we  see  that  the  numerator  in  the  expression  is 
negative  while  the  denominator  +  -0  . 

Therefore , 


x2  ♦  + 


+  OD 


(3.32) 


and  c  (real). 


Thus  stability  has  been  achieved  for  all  a  by  an  appropriate  choice 
of  Y  . 


3 . 5  Complex  Boundary  Condition  and  Boundary  Layer  Stability 

We  have  seen  that  the  stability  of  flows  In  which  the  velocity 
is  piecewise  constant  can  be  dramatically  affected  by  the  use  of  complex 
(in  Fourier  space)  boundary  conditions.  Since  our  main  interest  is  the 
stability  of  boundary  layer  flows,  in  which  the  velocity  varies 
continuously,  we  would  like  to  have  some  indication  of  what  can  bo  done 


We  consider  the  Rayleigh  equation  for  a  flow  U(z)  for  which 


U(0)  -  0  ,  U(*»)  ■  Uq  and  U'(z)  >  0  .  If  we  denote  the  stream 

function  by  $(z)eio^x~cC^  and  define 

♦  “  (J*)  *1 " e<iz*  h  m  ^ eOZ ♦'  (3,33) 

the  Rayleigh  equation  may  be  written 

t*-|-(i!)*sK7  (??)]♦ 

with  the  boundary  condition 


*<-)  “  (_J)  •  (3.35) 

The  boundary  condition  at  the  wall  (z  *  0)  will  be  specified  later. 
Since  U(z)  is  monotone,  and  has  a  finite  range,  It  makes  sense  to  use 
It  as  an  Independent  variable.  Then 

30-  *  ■  [it  (i  l)  +  onr-cT-  (i  o)J  *  (3,36> 

Now  U'  and  U' '  are  to  be  thought  of  as  functions  of  U  . 


We  of  course  cannot  solve  this  equation  In  general.  Let  us 
therefore  consider  the  long  wavelength  limit,  a  ♦  0  .  We  will  see  that 
It  Is  easy  to  arrange  the  wall  boundary  conditions  so  that 
ac  ♦  constant  In  this  limit.  The  Rayleigh  equation  now  reads 
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(3.37) 


1 


it  ■  (°  °L 

dU  ac  \l  Or 


and  its  solution  (incorporating  the  boundary  condition  at  z  ■  «•  )  is 


it'(U)  »  exp 


(-!) 

I  '  +  /  dv  ^irj 


(3. 38> 


The  boundary  condition  at  the  wall  is  expressed  in  general  as 


’*’2(0)  *'(0)  <T(  a,  ac) 


<|»  (0)  a<ji(0) 


(3.39) 


where  S  has  to  satisfy  certain  reality  and  causality  conditions.  A 
form  which  is  general  enough  for  our  purpose  is 


(3.40) 


where  a  and  b  are  dimensionless.  >  0  and  we  have  introduced 

D 

U'(0)  to  nondimensionalize  ac 


The  eigenvalue  condition  for  ac  is  now 


-I 


IT  '/U' 

dv  - 

ac 


1 


ac 

FToy 


(3.41) 


The  Integral  in  the  exponent  can  be  done  and  we  obtain  [recalling  that 
for  a  boundary  layer  flow  U'(")  ■  0  J, 


t 


il  1-9 


1 


(3.42) 


If,  for  simplicity,  we  set  a  -  0  ,  we  have 


ac  _  _ _ 1 

U' (0)  ~  in(l  +  lb)  ’ 


b  arbitrary. 


(3.43) 


The  imaginary  part  of  ac  may  be  made  as  large  and  negative  as  we  like 
by  taking  b  small  and  negative. 


What  does  this  mean?  A  boundary  layer  with  definite  sign  of 
U"  Is  Invlscldly  stable  anyway.  On  the  other  hand,  the  decay  rate  Is 
expected  to  go  to  zero  as  a  ♦  0  ,  indicating  a  nearby  Instability. 

With  the  new  boundary  condition  we  can  push  this  Instability  away.  A 
profile  with  an  inflection  point  may  be  inviscidly  unstable  for  small 
enough  a  .  We  apparently  can  control  this  instability  as  well.  We 
really  want  to  know  whether  we  can  improve  stablity  uniformly  for  all 
a  ,  and  the  above  discussion  gives  only  a  hint  that  this  is  possible. 


3.6  Conclusion 

From  these  simple  models  we  conclude  It  should  be  possible  by 
constructing  appropriate  boundary  materials  to  dramatically  affect 
stability. 
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BOUNDARY  LAYER  PERTURBATION  THEORY 


We  would  like  to  pose  a  problem  which  arises  in  the  study  of 
the  flat  plate  boundary  layer.  The  basic  equations  (the  Prandtl 
equations)  for  the  velocity  components  (  u  is  parallel  to  the  plate, 
v  ,  perpendicular)  are 


3u  3v 
3x  3y 


0 


3u 

3x 


+  V  -T— 


3u 

3y 


with  the  boundary  conditions 


u(0)  »  v(0)  *  0  ,  u(®)  =  Uq  v(»)  =  0 

This  system  is  parabolic:  given  an  initial  profile  u(xQ,y)  one 
integrates  forward  in  x  to  obtain  u(x,y)  at  any  later  x.  The 
standard  Blastus  solution  is  obtained  by  imposing  the  condition  of  self- 
similarity 


u(x,y)  »  F(y/u0/ \raj  , 

a  condition  which  can  be  met  because  of  the  Reynolds  scaling  property  of 
the  equations.  Suppose  the  initial  conditions  deviate  from  those 
required  to  produce  a  self-similar  solution.  What  happens  then?  Does 
the  solution  eventually  converge  to  Blasius,  or  does  it  perhaps  evolve 


away  from  it? 


To  attack  this  question  we  will  try  to  pose  the  linear 
stability  problem  for  small  perturbations  around  the  Blasius  solution  to 
the  Prandtl  equations.  The  equations  are  best  stated  in  terms  of  the 
stream  function  4)  , 


“-•$  '--U 


(4.1) 


A  self-similar  profile  must  result  from  a  stream  function  of  the  form 


-  AJTJx  fQ(n)  n-  y 


The  Prandtl  equations  then  reduce  to  an  ordinary  differential  equation 


f"'  +  2f  f"  -  0 
o  00 


(A. 2) 


f  (0)  -  f'(0)  -  0,  f'(«)  -  1 
00  o 


whose  solution  can  be  obtained  numerically. 


We  want  to  study  the  time  evolution  of  small  perturbations  on 
f.«  To  that  end  we  consider  stream  functions  of  the  general  form 


♦  -  /w*  (f0<n>  + 


(A. 3) 


and  expand  the  Prandtl  equations  to  first  order  In  fj  .  We  find 


I 


* 


t 


J 


I 


J 


I 


t 

* 


v  fr  -  »f;  n  *  c*  x)fo'f  i  - 0  <4-4) 

This  is  to  be  solved  subject  to  boundary  conditions  guaranteeing  that 
does  not  change  the  basic  flow: 

fj(0)  -  f'(0)  -  0,  f'(«0  -  0  (4,5) 


For  the  moment  the  parameter  a  ,  which  governs  the  downstream  growth 
rate  of  the  perturbation,  is  arbitrary.  In  fact,  the  differential 
equation  plus  boundary  conditions  pose  a  sort  of  eigenvalue  problem 
for  a  .  If  we  find  that  Re(a)  <  0  ,  then  we  will  conclude  that 
perturbations  on  the  standard  Blasius  profile  die  away  downstream,  a 
physically  satisfactory  result. 


To  see  whether  the  boundary  conditions  can  be  met,  we  look  at 
the  large  y  limit  of  Eq.  (4.4).  That,  in  turn,  is  governed  by  the 
large  y  limit  of  fQ(n)  •  T^e  asymptotic  behavior  of  fQ  is  found 
from  a  study  of  Eq.  (4.2)  to  be 


f  <n) 
o 


n-1.72  +  0(e~n2/4) 


(4.6) 


Substituting  this  into  Eq.  (4.4)  and  dropping  terms  of 
obtain 


we 


f'"  +  j  (n-1.72)  f''  -  af ^  -  0  . 
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This  is  a  standard  hyper geometric  equation  for  f1  with  two  non-trivial 
solutions.  The  three  possible  asymptotic  behaviors  of  fj  are 

2 

f  i  I  2)a+1/2  _JV  .  2. -(a+1)  ,,  7. 

fj  ~  1,  )  ,  e  (n  )  (*•?) 

All  solutions  satisfy  the  boundary  conditions  of  Eq.  (A. 5)  so  long  as 
a  +  1/2  <  1  .  (We  will  shortly  see  that  we  must  actually  impose  the 
condition  a  +  1/2  <  0)  .  For  a  less  than  this  limit,  any  a  is 
acceptable.  For  a  greater  than  this  limit,  only  discrete  values  of 
a  for  which  the  coefficient  of  the  (r^)a+^^  solution  vanishes  will  be 
acceptable. 


To  see  why  we  must  impose  the  condition  a  +  1/2  <  0  ,  it  is 
helpful  to  look  at  the  two  Important  definitions  of  boundary  layer 
thickness-displacement  thickness  (5^  and  momentum  thickness  (5^)  • 
The  definitions  are 


62  “  f  dy  TT  (l  "  T") 

O  «  « 


(A. 6) 


where  u  is  the  component  of  velocity  parallel  to  the  wall.  6^  is 
equal  to  the  outward  displacement  of  the  stream  lines  far  from  the  wall 
compared  to  their  position  when  the  wall  is  absent.  62  is  the 
thickness  which  would  have  to  be  removed  from  the  free  stream  flow  in 
order  to  match  the  momentum  actually  lost  to  the  wall.  If  we  use  our 
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-  -J.  ' 


general  expression  for  <jj  [Eq.  (A. 3)],  remember  that  u  ■  3<|//3x  , 
expand  to  first  order  In  fj  and  Integrate  by  parts  whenever  necessary, 
using  the  equations  satisfied  by  fQ  and  fj  ,  we  obtain 

“  [1*72  +  xa(-fj(-))] 
f ' ' (0) 

52  -  .66  +  xa  a|l/2  (4*7) 

Since  the  displacement  thickness  must  be  finite  on  physical  grounds,  it 
is  apparent  that  we  must  impose  the  condition  f^(»>)  <  *  .  By  the 
previous  paragraph,  this  means  that  we  must  choose  a  +  1/2  <0  or  find 
a  (hypothetical)  discrete  eigenvalue  of  a  for  which  the  offending 

(n^)a+^^  term  does  not  appear  in  the  asymptotic  behavior  of  f^  • 
Since  a  numerical  search  did  not  show  any  sign  of  the  existence  of  a 
(real)  eigenvalue,  we  shall  henceforth  ignore  this  possibility. 


t 


» 


As  has  been  noted  in  many  places,  the  parameter  which  is  of  the 
greatest  importance  for  stability  is  the  ratio  6^/52  •  A  10%  change  in 
this  parameter  can  cause  an  order  of  magnitude  change  in  the  critical 
Reynolds  number.  To  first  order  in  f^  we  have 


a 


x 


c( 


#(a) 

1777 


*(a) 


fi(“) 

f''(0) 


c  -  f('(0) 


» 
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f  j ' (0)  can  be  thought  of  as  the  constant  normalizing  fjj  one 
Integrates  Eq.  (4.4)  with  the  starting  conditions 
fj(0)  -  fj(0)  -  0  f''(0)  ■  c  .  The  quantity  *(a)  is  obtained  by 
Integrating  the  differential  equation  out  to  Infinity,  and  is  Just  a 
function  of  a  ,  defined  for  a  <  -1/2  . 

We  have  done  the  requisite  numerical  integration  for  a  range  of 
values  of  a  <  -1/2  and  we  find  that  by  choosing  the  sign  of  c 
properly  we  can  always  arrange  the  perturbation  to  decrease  the  quantity 
6^2  thereby  Increasing  stability.  Because  a  <  -1/2  ,  the 
perturbation  must  die  away  downstream  at  least  as  rapidly  as  x*1/2  . 
Since  X-*/2  is  a  rather  slow  rate  of  decrease,  there  might  be  some 
advantage  to  attempting  to  set  up  the  perturbation  for  which  a  ■  -1/2 
— the  stability  Increase  which  Is  achieved  then  decays  very  slowly, 
perhaps  slowly  enough  to  maintain  laminar  flow  over  a  region  of 
Interesting  size. 

The  issue  we  have  not  yet  explored  Is  the  question  of  whether 
there  may  be  complex  discrete  eigenvalues  and  whether  to  represent  a 
general  perturbation,  complex  values  of  a  are  required.  We  would  like 
to  come  back  to  this  question  In  the  future,  although  we  expect  It  to  be 
a  rather  difficult  question  to  analyze,  since  the  operator  Involved  is 
not  hermitean. 
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APPENDIX 


Here  we  want  to  show  that  there  Is  an  unstable  node  for  a 

slightly  smaller  than  a  .  It  was  shown  in  Eq.  (2.23)  that 

s 


/ 


2  Z1 


..2  +  U"  m2L, 

>  I  +XJ^  jdz  * 


U(z2)  I 
~t: — 


I  <J»(Z  x )  I  ‘ 
— 5T~ 


c 


I  <J>|  2d  z 


(Al) 


Consider  the  change  of  -  a  with  c  In  the  vicinity  of  cs  •  Since 
2 

-  Is  stationary  we  merely  need  to  differentiate  with  respect  to  the 

explicit  occurrence  of  c  .  Since  we  are  looking  for  unstable  modes  we 

evaluate  at  c  «  c  +  ie 
s 


3a 

3c 


f  2  U"  l*l' 

ix  (U-c-ie)‘ 


dz 


r*  : 

J  l+l 


(A2) 


dz 


Since  U"(z_)  *  0  the  integrand  in  the  numerator  of  Eq.  (A2)  has  a 

d 

simple  pole  at  z  ■  zg  .  The  ic  tells  how  we  are  to  go  around  this. 
We  find 


A- 1 


*  <«- 


f*2  u"i*i2dj  ,  ’‘l**"  u'(«.)K"(s)|^.)|' 

J  [B.(|J|, - 


I ♦(z ) I  dz 


1*6*9 


-2a  wl  “  a  +  lb  • 

°8 


therefore, 


3c l  -2au(a-ib) 

3^1  ■  — y-  2 

a  a  +  b 

8 


/2  2 
U"  !  ♦!  dz 

z  U  -  U(z  ) 
1  s 


I ♦Cz ) |  dz 


TrU* '  )| 2  sgn  U'(zs) 

(u'(zs)]2  /  I ♦(*) I 2dz 


Then  for  small  a  -  a 


A- 2 


2a  b 
s 

"2  .2 
a  +  b 


(a  ~  afl )  . 


Let  us  see  how  this  checks  with  our  example. 

U(z)  »  sin  z,  z  -  0  . 

s 

Then  U'(z  )  -  +  1  ,  U"'(z  )  -  -1  . 

s  s 

Therefore,  b  <  0  . 

^  a  >  ag  E<1*  (A®)  gives  c^  <  o  ,  contrary  to  our  assumption.  On  the 
other  hand  if  a  <  ag  we  see  we  do  indeed  obtain  a  >  0  .  We 
conclude  there  is  an  unstable  mode  for  a  slightly  smaller  than  ag  ! 
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